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Abstract 

We consider the classical evolution of a particle on a graph by using 
a time-continuous Frobenius-Perron operator which generalizes previous 
propositions. In this way, the relaxation rates as well as the chaotic prop- 
erties can be defined for the time-continuous classical dynamics on graphs. 
These properties are given as the zeros of some periodic-orbit zeta func- 
tions. We consider in detail the case of infinite periodic graphs where the 
particle undergoes a diffusion process. The infinite spatial extension is 
taken into account by Fourier transforms which decompose the observables 
and probability densities into sectors corresponding to different values of 
the wave number. The hydrodynamic modes of diffusion are studied by 
an eigenvalue problem of a Frobenius-Perron operator corresponding to a 
given sector. The diffusion coefficient is obtained from the hydrodynamic 
modes of diffusion and has the Green-Kubo form. Moreover, we study 
finite but large open graphs which converge to the infinite periodic graph 
when their size goes to infinity. The lifetime of the particle on the open 
graph is shown to correspond to the lifetime of a system which undergoes 
a diffusion process before it escapes. 

PACS numbers: 

02.50.-r (Probability theory, stochastic processes, and statistics); 
03.65.Sq (Semiclassical theories and applications); 
05.60.Cd (Classical transport); 

45.05.-l-x (General theory of classical mechanics of discrete systems). 



1 Introduction 

The study of classical dynamics on graphs is motivated by the recent discovery 
that quantum graphs have similar spectral statistics of energy levels as the 
classically chaotic quantum systems ^ . Since this pioneering work by Kottos 
and Smilansky, several studies have been devoted to the spectral properties of 
quantum graphs |^. ||]. However, the classical dynamics, which is of great 
importance for the understanding of the short-wavelength quantum properties. 
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has not yet been considered in detail. In Refs. iQ, ||], a classical dynamics has 
been considered in which the particles are supposed to move on the graph with 
a discrete and isochronous (topological) time, ignoring the different lengths of 
the bonds composing the graph. 

The purpose of the present paper is to develop the theory of the classical 
dynamics on graphs by considering the motion of particles in real time. This 
generalization is important if we want to compare the classical and quantum 
quantities, especially, with regard to the time-dependent properties in open or 
spatially extended graphs. A real-time classical dynamics on graphs should 
allow us to define kinetic and transport properties such as the classical escape 
rates and the diffusion coefficients, as well as the characteristic quantities of 
chaos such as the Kolmogorov- Sinai (KS) and the topological entropies per unit 
time. 

An important question concerns the nature of the classical dynamics on a 
graph. A graph is a network of bonds on which the classical particle has a one- 
dimensional uniform motion at constant energy. The bonds are interconnected 
at vertices where several bonds meet. The number of bonds connected with a 
vertex is called the valence of the vertex. A quantum mechanics has been defined 
on such graphs by considering a wavefunction extending on all the bonds ^ . 
This wavefunction has been supposed to obey the one-dimensional Schrodinger 
equation on each bond. The Schrodinger equation is supplemented by boundary 
conditions at the vertices. The boundary conditions at a vertex determine the 
quantum amplitudes of the outgoing waves in terms of the amplitudes of the 
ingoing waves and, thus, the transmission and reflection amplitudes of that 
particular vertex. 

In the classical limit, the Schrodinger equation leads to Hamilton's classical 
equations for the one-dimensional motion of a particle on each bond. When a 
vertex is reached, the square moduli of the quantum amplitudes give the proba- 
bilities that the particle be reflected back to the ingoing bond or be transmitted 
to one of the other bonds connected with the vertex. In the classical limit of 
arbitrarily short wavelengths, the transmission and reflection probabilities do 
not reduce to the trivial ones (i.e., to zero and one) for typical graphs. Ac- 
cordingly, the limiting classical dynamics on graphs is in general a combination 
of the uniform motion of the particle on the bonds with random transitions at 
the vertices. This dynamical randomness which naturally appears in the clas- 
sical limit is at the origin of a splitting of the classical trajectory into a tree of 
trajectories. This feature is not new and has already been observed in several 
processes such as the ray splitting in billiards divided by a potential step 
or the scattering on a wedge 0. We should emphasize that this dynamical 
randomness manifests itself only on subsets with a dimension lower than the 
phase space dimension and not in the bulk of phase space, so that the classical 
graphs share many properties of the deterministic chaotic systems, as we shall 
see below. 

The dynamical randomness of the classical dynamics on graphs requires a 
Liouvillian approach to describe the time evolution of the probability density 
to find the particle somewhere on the graph. Accordingly, one of our first goals 
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below will be to derive the Frobenius-Perron operator as well as the associated 
master equation for the graphs. This operator is introduced by noticing that 
the classical dynamics on a graph is equivalent to a random suspended flow 
determined by the lengths of the bonds, the velocity of the particle, and the 
transition probabilities. 

A consequence of the dynamical randomness is the relaxation of the prob- 
ability density toward the equilibrium density in typical closed graphs, or to 
zero in open graphs or in graphs of infinite extension. This relaxation can be 
characterized by the decay rates which are given by solving the eigenvalue prob- 
lem of the Frobenius-Perron operator. The characteristic determinant of the 
Frobenius-Perron operator defines a classical zeta function and its zeros - also 
called the PoUicott-Ruelle resonances - give the decay rates. The leading decay 
rate is the so-called escape rate. The Pollicott-Ruelle resonances have a par- 
ticularly important role to play because they control the decay or relaxation 
and they also manifest themselves in the quantum scattering properties of open 
systems, as revealed by a recent experiment by Sridhar and coworkers The 
decay rates are time-dependent properties so that they require to consider the 
continuous-time classical dynamics to be defined. 

Besides, we define a continuous-time "topological pressure" function from 
which the different chaotic properties of the classical dynamics on graphs can 
be deduced. This function allows us to define the KS and topological entropies 
per unit time, as well as an effective positive Lyapunov exponent for the graph. 

We shall also show how diffusion can be studied on spatially periodic graphs 
thanks to our Frobenius-Perron operator and its decay rates. Here, we consider 
graphs that are constructed by the repetition of a unit cell. When the cell is 
repeated an infinite number of times we form a periodic graph. Such spatially 
extended periodic systems are interesting for the study of transport properties. 
In fact at the classical level it has been shown in several works that relationships 
exist between the chaotic dynamics and the normal transport properties such 
as diffusion |^ and the thermal conductivity [0, which have been studied in 
the periodic Lorentz gas. In the present paper, we obtain the continuous-time 
diffusion properties for the spatially periodic graphs. Moreover, we also study 
the escape rate in large but finite open graphs and we show that this rate is 
related, on the one hand, to the diffusion coefficient and, on the other hand, to 
the effective Lyapunov exponent and the KS entropy per unit time. 

The plan of the paper is the following. Sec. || contains a ge neral introduction 



to the graphs and their classical dynamics. In Subsec. 2.2, we introduce the 
evolution using the aforementioned random suspended flow and, therefore, we 
can follow the approach developed in Ref. for the study of relaxation and 
chaotic properties at the level of the Liouvillian dynamics which is developed 



in Sec. H. The Frobenius-Perron operator is derived in Subsec. 3.1, In Subsec. 



3.2, we present an alternative derivation of the Frobenius-Perron operator and 
its eigenvalues and eigenstates, which is based on a master-equation approach, 
familiar in the context of stochastic processes. Both approaches are shown 
to be equivalent. In Sec. ^, we study the relaxation and ergodic properties 
of the graphs in terms of the classical zcta function and its Pollicott-Ruelle 
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resonances. In Sec. the large-deviation formalism is introduced which allows 
us to characterize the chaotic properties of these systems. In Sec. ^, the theory 
is applied to classical scattering on open graphs. The case of infinite periodic 
graphs is considered in Sec. ^, where we obtain the diffusion coefficient and we 
show that it can be written in the form of a Green-Kubo formula. In Sec. ||, 
we consider finite open graphs of the scattering type, where the particle escape 
to infinity, and we show how the diffusion coefficient can be related to the 
escape rate and the chaotic properties. The case of infinite disordered graphs is 
considered in Sec. 0. Conclusions are drawn in Sec. O. 



2 The graphs and their classical dynamics 
2.1 Definition of the graphs 

Now, let us introduce graphs as geometrical objects where a particle moves. 
Graphs are V vertices connected by B bonds. Each bond b connects two vertices, 
i and j. We can assign an orientation to each bond and define "oriented or 
directed bonds". Here one fixes the direction of the bond [i.j] and call b = 
the bond oriented from itoj. The same bond but oriented from j to i is denoted 

b = {j, i). We notice that b — b. A graph with B bonds has 2B directed bonds. 
The valence Vi of a vertex is the number of bonds that meet at the vertex i. 

Metric information is introduced by assigning a length lb to each bond b. In 
order to define the position of a particle on the graph, we introduce a coordinate 
Xb on each bond b = We can assign either the coordinate or 

The ffist one is defined such that a;(ij) = at i and = lb at j, whereas 

^U,i) = at J and xq- ^) = at i. Once the orientation is given, the position of 
a particle on the graph is determined by the coordinate Xb where < a;;, < 4- 
The index b identifies the bond and the value of Xb the position on this bond. 

For some purposes, it is convenient to consider b and b as different bonds 
within the formalism. Of course, the physical quantities defined on each of 
them must satisfy some consistency relations. In particular, we should have 
that = lb and = h — Xb- 

A particle on a graph moves freely as long as it is on a bond. The vertices are 
singular points, and it is not possible to write down the analogue of Newton's 
equations at the vertices. Instead we have to introduce transition probabilities 
from bond to bond. These transition probabilities introduce a dynamical ran- 
domness which is coming from the quantum dynamics in the classical limit. In 
this sense, the classical dynamics on graphs turns out to be intrinsically random. 

The reflection and transmission (transition) probabilities are determined by 
the quantum dynamics on the graph. This latter introduces the probability 
amplitudes Tbb' for a transition from the bond b' to the bond b. We shall 
show in a separate paper |pT| that the random classical dynamics defined in 
the present paper, with the transition probabilities defined by Pbb' = \Tbb'\^ is, 
indeed, the classical limit of the quantum dynamics on graphs. For example. 



4 



we may consider a quantum graph with transition amphtudcs of the form 

Tbt' = Cbt' ( -<5b.6- ) (1) 

where Cbb' is 1 if the bond b' is connected with the bond 6 and zero otherwise 
and i>bb' is the valence of the vertex that connects b' with b. Such probability 
amplitudes are obtained once we impose the continuity of the wave function and 
the current conservation at each vertex. 

For the classical dynamics on graphs, the energy of the particle is conserved 
during the free motion in the bonds and also in the transition to other bonds. 
Accordingly, we consider the surface of constant energy in the phase space de- 
termined by the coordinate of the particle, that is which specifies a bond and 
the position with respect to a vertex. The momentum is simply given by the 
direction in which the particle moves on the bond. The modulus of the momen- 
tum is fixed by the energy. We see that position and direction can be combined 
together if we speak about position in a given directed bond. In this way our 
phase space are the points of the directed bonds. The equation of motion is 
thus 



— = i; = ^/2E/M , for < x = Xb < h , (2) 

where v is the velocity in absolute value, E the energy, and M the mass of the 
particle. When the particle reaches the end Xf,' — lb' of the bond b' a transition 
can bring it at the beginning Xb = oi the bond b. According to the above 
discussion, we assume moreover that this transition from the bond b' to the 
bond b has the probability Pbb' to occur: 

transition b' ^ b with probability Pbb' (3) 

By the conservation of the total probability, the transition probabilities must 
satisfy 

T.P"'^'-^ (4) 

b 

which means that the vector {1,1,... ,l}is always a left eigenvector with eigen- 
value 1 for the transition matrix P = {Pbb'}b%=i- 

We may assume that the system has the property of microscopic reversibility 
(i.e., of detailed balancing) according to which the probability of the transition 
6' ^ 6 is equal to the probability of the time-reversed transition b ^ b': Pbb' = 
P^-,g, as expected for instance in absence of a magnetic field. As a consequence 
of detailed balancing, the matrix P is a bi-stochastic matrix, i.e., it satisfies 

Pbb' = X^b' ^bb' = 1, whereupon the vector {1,1,... , 1} is both a right and 
left eigenvector of P with eigenvalue 1. This is the case for a finite graph with 
transition probabilities Pbb' = iTfcfe'P given by the amplitudes (|^). 
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2.2 The classical dynamics on graphs as a random sus- 
pended flow 

The description given above is analogous to the dynamics of a so-called sus- 
pended flow In fact, we can consider the set of points {xb = 0, V6}, i.e., the 
set of all vertices, as a surface of section. We attach to each of these points a seg- 
ment (here the directed bond) characterized by a coordinate < x < lb- When 
the trajectory reaches the point x — lb it performs another passage through the 
surface. Thus the flow is suspended over the Poincare surface of section made 
of the vertices in the phase space of the directed bonds. 

For convenience, instead of the previous notation Xb, the position (in phase 
space) will be referred to as the pair [b, x] where b indicates the directed bond 
and X is the position on this bond, i.e., < x < lb- 

A realization of the random process on the graph (i.e., a trajectory) can be 
identified with the sequence of traversed bonds • • • 6-2&-if'o^i&2 • • ■ (which is 
enough to determine the evolution on the surface of section). The probability 
of such a trajectory is given by • • • Pb2biPbiboPbob^iPb_ib_2 ■ ■■ ■ 

An initial condition [feoi^:] of this trajectory is denoted by the dotted bi- 
infinite sequence • • • 6_2&-i • 60^1^2 • • • together with the position < x < Ibg. 
For a given trajectory, we divide the time axis in intervals of duration 
extending from 

_j_ _p . . . ^ _j_ ^ . . . ^ -|- 

V V V V V V V V V 

where v is the velocity of the particle which travels freely in the bonds. At each 
vertex, the particle changes its direction but keeps constant its kinetic energy. 

For a trajectory p that, at time i = 0, is at the position [bo, x] we define the 
forward evolution operator with i > by 

^p[bQ,x] ^ [bo,vt + x] if 0<x + vt<ho (5) 

i.e., the evolution is the one of a free particle as long as the particle stays in the 
bond 6o, and 

[5o, x] = [5„, x + vt~ h„-i - 4„_2 ho] 

if < X + vt ~ h„^i - h„-2 — ■ ■ ■ - ho < h„ (6) 

which follows from the fact that, for the given trajectory p, the bond and the 
position where the particle stands at a given time is fixed by the lengths tra- 
versed at previous times and by the constant velocity v. Analogously, we also 
introduce a backward evolution for t < Q: 

$p 1*1 [^0,2^] = [&o,a; - wlil] if Q<x-v\t\<ho 

and 

$-1*1 [bo, x] = [6_„, X - v\t\ + h^„ + 4_„+i + ■ • • + h.i] 

if < x - v\t\ + h-n + 4_„+i H h /fc_i < /b_„ 
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3 The Liouvillian description 



3.1 The Frobenius-Perron operator 

On the graph, we want to study the time evolution of the probabihty density 
p ( [6, x\ ,t). This density determines the probabihty p ( [6, x\ , t) dx of finding the 
particle in the bond h with position in [x, x + dx\ at time t. 

For a general Markov process, the time evolution of the probability density 
is ruled by the Chapman-Kolmogorov equation 

p(C,t) = j diaV {iAioM P{^M (7) 

where V tl^Oj ^o) is the conditional probability density that the particle be in 
the state ^ at time t given that it was in the state at the initial time . This 
conditional probability density defines the integral kernel of the time-evolution 
operator, which is linear. The conditional probability density can be expressed 
as a sum (or integral) over all the paths joining the initial state to the final one 
within the given lapse of time. 

In the case of graphs, each state is given by an directed bond and a position 
on this bond: ^ = [6, x]. A path or trajectory is a bi- infinite sequence of 
directed bonds as described in the preceding section. As soon as the path or 
trajectory p is known, the sequence of visited bonds is fixed so that the motion 
is determined to be the time translation at velocity v given by Eqs. d)-®. In 
this case, the conditional probability density of finding the particle in position 
[6, a;] at time t given it was in [5o,a:o] at the initial time io = is given by a 
kind of Dirac delta density 5 ([6, x] — <&* [6o, xq\) . Along the path p, the particle 
meets several successive vertices where the conditional probability density is 
expressed in terms of the conditional probability to reach the final bond h = hn 
within the time given the initial condition [bp, xq\: 

We notice that the integer n is fixed by the trajectory p, the initial condition 
[&o, 2^0], and the elapsed time t. The number of the path probabilities (H) which 
are non- vanishing is always finite if there is no sequence of lengths accumulating 
to zero for the graph under study. 

For a graph, the kernel of the evolution operator is thus given by 

V{[b, x] , t\ [&o, xo] , 0) = ^ Pp{t, [60, xo]) S {[b, x] - [bo, xo]) , (9) 
{p} 

where the sum is performed over a finite number of paths. By analogy with Eq. 
(0) the density is given by 

p{[b,x],t) ^J2f I JZ-^P*^^' [^0,2:0]) S {[b,x] - $* [bo,xn]) ) p ([feo, xq] , 0) 

bo \{p} 
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and integrating the Dirac delta densit}^ we finally get 

{p} 

where the sum is over all the trajectories that go backward in time from the 
current point [b, x] . 

In this way, we have defined the Frobenius-Perron operator P* as 

F [b, x]^j2 ^]) ^ ('^p' ^0 (11) 

{p} 

where {F[b, x]} is a vector of 2B functions defined on the directed bonds. 

We now turn to the determination of the spectrum of the Frobenius-Perron 
operator. With this aim, we take the Laplace transform of the Frobenius-Perron 
operator given by Eq. (O): 



/■oo poo 

/ e-'^P* F[b,x]dt^ dte-^*^Fp(t,$p' P($p* 
^0 Jo 

In order to evaluate the Laplace transform, we have to decompose the sum over 
the paths {p} into the different classes of terms corresponding to paths in which 
n = 0, 1, 2, ... bonds are visited during the time t. This decomposition leads to 

V f e-"'Pp F [b, x]) dt^ /" " e-'*F [6, x - vt] dt 

n 



{p} 



i=l 



dt 



with Pp-n — Pbb-iPb-ib^2 ' ' ' Pb-n.+ib-n the probability of a path and 
^lp_„j the sum over these trajectories. A change of variable transforms the 
previous equation in the following 



V/ e-"*PpF($;*[5,2;]) = ie-"" ( / e'"^ F [b,x']dx' 

+ E E n^^^+i"-. / <i-'^F[b.n.Adx') (12) 

n=l{p-"}V»=l /-^O J 

where we introduced the quantity 

Qbb'{s)^Pbb' G-''-^ (13) 



^ We emphasize the analogy with deterministic processes for which 'P{x,t\xQ,0) 
&[x- (/>*(xo)] so that p{x, t) = I I p [4>-\x), 0] . 
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and here we identify bo with b. Now we have to perform the sum over all the 
realizations in the right-hand side of Eq. This is a sum over all the 

trajectories 6_i6_2 ■ ■ -b-n+ib^n which leads to the formation^ of the matrix 
Q", i.e., Q raised to the power n. Accordingly, we get 

f e-''*P* F [b, x]dt^- c-'''^'" ( / e'^F [b, x'] dx' 
Jo ^ Uo 

+ E E (^")bb „ / e^"^^ dx' \ . (14) 

n=lfc-,. J 

We define the vector f (s) = {/i(s), /2b(s)} with the components of this 
vector given by the functions 

/b(s)= [\'"^F[b,x']dx' (15) 
Jo 

The matrix Q(s) acts on these vectors f(s) through the relation 

iQ-f\i.s) = J2Qwis)fb'is). (16) 

b' 

This matrix can be interpreted as the Frobenius-Perron operator of the evo- 
lution reduced to the surface of section Q . This matrix depends on the Laplace 
variable s which will give the relaxation rate of the system. With these defini- 
tions, we can write 

EE(Q"W,. / e^vF[6_„,x']d^' = E(Q"-f)6(^)= r^-f 



where we used the relation X^^^i Q" = ^ consequence, Eq. ( p4[ ) 

becomes 

(17) 

We are now at a few steps from determining the eigenvalues (and eigenvec- 
tors) of P*. This is done by first studying the solutions of 

Q{s) ■ f (,s) - {(s). 

These solutions exist only if s belongs to the (complex) set {sj} of solutions of 
the following characteristic determinant 

dot [I - Q(s)] = (18) 

'^T^b_i---b_„ Qbob^i ■ ■ ■Q6_„+ib_„ = Ei,_„(Q")fco6-Ti 
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We denote these particular vectors by Xj ^^^d their components by Xj [6] with 
6=1, 2B, whereupon 

Qi-s,)-Xj=Xr (19) 
The left vector which is adjoint to the right vector Xj is given by 

Qis,y-Xj=Xj- (20) 

The relation to the eigenvalue problem of the flow is established as follows. 
Suppose that 'i'j [b, x] is an eigenstate of P* with eigenvalue e^^ * [for the moment 
Sj is not determined but we call it this way because it will turn out to be one 
of the solutions of Eq. (|l|)], i.e., 

P**j = e''^**j [6,a;] (21) 

with Re Sj < because the density is not expected to increase with the time. 
For the forward semigroup, the zeros of Eq. ( JTsj ) are expected in the region 
Re Sj < 0. 



Taking the Laplace transform of Eq. (|2l|), we get 

/•oo poo 

/ dt e""* P* -fj [b, x]^ dt e'^'^-'^'-fj [b, x] = 
Jo Jo 



S — So 



■=3 

If we introduce the vector T(s) — {Ti(s), T2b(s)} defined by the components 

Tf,(s) = f dxe'^^j [b,x] 
Jo 

and use the same calculation that led to Eq. (|l7|), the eigenvalue equation 
becomes 

Q _ V e'^ [b,x] 



dx'e^^^,[b,x']+{-^-r] {s) = (22) 



for < X < lb- Setting a; = in Eq. (^2|), we have 

{s~s,)Q{s)-T{s)=v [l-Q(s)].*, 
with the vector = {^j[6, 0]}^;f^. For s = Sj, we get that 

which shows that sj is a solution of Eq. as we anticipated and that the 

eigenstate of the flow at x ~ 0, [b, 0], may be identified with the vector which 
is solution of Eq. (|l9|): 

^,[b,0]^Xj[b]^J2Qbb'{s,)xAb'] 
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To determine the eigenstates of the flow for the other values of x we differentiate 
Eq. ( p2| ) with respect to x and we get 

e^^ ^, %x] = i£!iliM + -^'-d'^^.[bM 



from which we obtain 



the integration of which gives 

*j [5, x] = e'^^-'^j [5, 0] = e-"^ " xj [b] for Q<x<lb. (23) 

The eigenstate increases exponentially along each directed bond. This expo- 
nential increase does not constitute a problem because the time evolution gen- 
erates the overall exponential decay of Eq. (|T]). Therefore, we see that the 
vectors Xj W\ which are solutions of Eq. (19) determine the eigenstates of the 
Frobenius-Perron operator. These eigenstates are very important in nonequilib- 
rium statistical mechanics since they provide the link between the microscopic 
and the phenomenological description of the system 

3.2 A master-equation approach 

In this section, we develop an alternative derivation of the results of the previous 
subsection but here by using a master equation. 

If at a given time t the particle is at the end of a bond, say [6',/fc'], the 
particle has to go instantaneously to another directed bond with probability 
Pbb', i.e., 

p{[h,Q],t) = Y.Pbb'P{[b'M.t) (24) 

b' 

Now, since the evolution is deterministic along the bonds [see Eq. (||)] we have 
that^ 



=p([&,0],t) 



and also 

pi\h'M.t) = p{[b'Jb' -vt'],t-t') 

Choosing t' — , we have 



This is because p{[b,x\ ,* + f ) = p{[b,x\ ,t) = p (4>"^ [b,x\ ,t] = p([b,0] ,t). 
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where x' is arbitrary. Hence, Eq. ( |24| ) becomes 

P ^] ' * + ~) ^^^Pwp[ [6', a:^6'] , 



Xh' 



where xv may be chosen arbitrarily on each bond 6'. With the replacement 
t + - ^ i, we finally obtain 



P 



(25) 



This is the master equation which rules the time evolution on the graph. It is 
a Markovian equation with a time delay. The master equation ( p5| ) differs from 
Eq. ( p^ ) in the sense that Eq. ( p5| ) relates the probability densities before and 
after the transitions although Eq. ( p^ ) relates the density at time t to the initial 
density through a varying number of transitions, depending on the path -p. 

Stationary solutions of the equation (p5|), satisfying /9([&, a;] ,i) — p([6, x]) 
for all t, exist if the matrix P as an eigenvalue 1. This is the case for closed 
graphs. For open graphs, the density decays in time in a way that we shall 
determine below. 

The master equation ( p5|) can be iterated. For instance the second iteration 
gives 



h'h" P 



b'b" 



b ,Xb', 



X + lb' + lb 



Xb'b" 



and in general 

p{[b,x],t)^ ^ Pbb'Pb'b" ■■■ Pb('^-i>b("> 

+ SLi 4<') - ^b'-bi'^) 



b'b"---U") 



p {y^\xb'...u^) 



(26) 



There exists an integer n for which we find (at least one) solution of 

X + lb" H h lb('^) - Xb'...b(^) 



with < Xb,...bM < Ib(^) (27) 



for some path b'b" ■ ■ ■ fe'"^ Accordingly, we split the sum ( [2^ ) in two terms, the 
first one with all the possible paths for which there exists a value Xbi...b{n) which 
solves Eq. (^^ with the smallest integer n (^' denotes the sum over these 
paths) and the other term containing the rest of (|26[), i.e.. 



p{[b,x\,t)^ ^ Pbb'Pb'b" ■■■Pb(^-^)b(^)p(^ 



,Xb,...b(^) 



,0 



+ 5Z ^bb'Pl 



b'b"---b<-"> 

'b'b" ■ ■ ■ Pb('^-i)b("'iP 



b'b"---b^'^) 



, Xbi...b(ri) 



^+Er=i4<') - Xb'-b('-)' 
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and we proceed iteratively with the second term, that is, we look for the smallest 
n for which there exists a path for which a solution of Eq. ( p7| ) exists and so 
on. Thus we finally have 



n b'b". ..;,(") 

with 

n 

a;^'. ..(,(") = .T - ft + ^ ■ 



b ,xy...^(„ 



,0) (28) 



Accordingly, p ( [6, x\ , t) is given by a sum over the initial conditions ^ , Xf,/ 



and over all the paths that connect 



, Xfj, 



with [b, x] in a time t. Each 



given path contributes to this sum by its probability multiplied by the proba- 



bility density p ^ 



b ,Xh,...h(,„) 



) . Using the notation introduced before [see 



Eq. (g)] Eq. ( Pq ) can be written as 

p([6,a;],t) = ^Pp(t,$;*[6,x]) p( 
{p} 



b ,Xi,'...b(^) 



,0) 



(29) 



with a;5/...f,c„) = a; — + X^iLi ^b(»)- We see that this equation coincides with 
Eq. (p^), which shows that both approaches are equivalent. In fact, if we 
write the sequence bb' ■ ■ ■ in the form 66_i • • • 6_„ and if we remember that 

$-*[6,a:] = [&_„, x - + ELi if < ^ - + YJl=i^b^^ < the 
equivalence is established. 

Now we turn to the determination of the spectrum of the Frobenius-Perron 
operator from the master equation. Taking the Laplace transform of Eq. ( p5| ) 
or equivalently of 

lb' + X 



we have 



\[bA,t-^)=Y.Pwp([b',0], 



^ e-^*'p ([6, 0] , t') dt' = Pbb' e-'^ J ^ e-'''p ([&', 0] , t') dt' (30) 



after some simple changes of variable.^ Considering the definition of Eq. (^ 
and that p ([6, 0] , t < 0) = 0, Eq. (|0) reads 



POO POO 

/ e-'*'p{[b,0],t')dt' =y2Qbb'{s) e-'''pi[b\0],t')dt'. (31) 
Jo ^, Jo 



At the left-hand side t' = t — x/v and at the right-hand side t' = t — x/v — ly /v. 
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oc 

St , 



Defining 

Jo 

Eq. becomes 

b' 

which has sohitions only if s belongs to the set {sj} of solutions of 

det [I - Q(s)] = 

and 

X.-Q(5.)-X, • (32) 

The eigenstates of the flow 

p,{[b,x],t)^e'''pj{[b,x],G) (33) 



are determined as follows. We replace Eq. in the master equation (g5|) 

from where we directly get 

{e^^^p, ,0)} ^Y.Qbb' {sj) {e^^i^p, ([6', a:] ,0)} . (34) 



Comparing Eq. (|3^) and Eq. (|3j) we have that the eigenstates of the Frobenius- 
Perron operator of the flow are given by 

Pj {[b, , 0) = G^^'^Xj W\ for < a: < 4 

and we have recovered the same results as previously obtained with the suspended- 
flow approach. 



4 The relaxation and ergodic properties 

4.1 The spectral decomposition of the Probenius-Perron 
operator 

Thanks to the knowledge of Frobenius-Perron operator, we can study the time 
evolution of the statistical averages of the physical observables A[6, x] defined 
on the bonds of the graphs as 

2S , 

{A)t^y^Y A[b,x]p{[b,xlt)dx={A\P'po) (35) 
b=i -^0 
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where po denotes the initial probabihty density and where we have introduced 
the inner product 



2-B , .1^ 

{F\G) = V 7- / F[b,x]* G[b,x]dx (36) 
&=i -^0 

between two vectors of 2B functions x] and G'[6, x] defined on the bonds. 
We have here used the fact that the Frobenius-Perron operator rules the time 
evolution of all the statistical averages. If the observable is equal to the unity, 
A = 1, the conservation of the total probability imposes the normalization 
condition = 1, which is satisfied by the Frobenius-Perron operator. 

If we are interested in the time evolution at long times and especially in the 
relaxation, we may consider an asymptotic expansion valid for t +00 of the 
form 

{A)t = {A\P'po) = Y.{A\^,) e^^-* (*,|po) + • • • (37) 



as a sum of exponential functions, together with possible extra terms such as 
powers of the time multiplied by exponentials i™exp(sjt). In this spectral de- 
composition, we have introduced the right and left eigenstates of the Frobenius- 
Perron operator 

= e'^^^j (38) 
P*t|-^. = e*l**j- (39) 

Since the Frobenius-Perron operator is not a unitary operator we should expect 
Jordan-block structures and associated root states different from the eigenstates. 
Such Jordan-block structures are known to generate time dependences of the 
form t"^ exp{sjt). We shall argue below that such time behavior is not typical 
in classical graphs. 

With the aim of determining the spectral decomposition (|37|), we take its 
Laplace transform: 

f°° 1 

/ e-^*(AMt = ^(A|vl/,-) + 

Jo J Sj 

which allows us to identify the relaxation rates —Sj with the poles of the Laplace 
transform and the eigenstates from the residues of these poles. For this purpose, 
we use the Laplace transform of the Frobenius-Perron operator given by Eq. 
(pi]), which we integrate with the observable quantity A [6, a:]. We get 

/ Q-"* {A\P^ pa) dt = V-p / dx dx'e''^A[h,x\pQ[b,x'] 
Jo J Jo Jq 
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where we introduced the vector a(s) of components 



ab(s) = — / e'"' A[b,x]dx (41) 
vlb Jo 

and where we used the definition ( p^ with the initial probabihty density F = pQ. 

In Eq. ( pO| ) , the first term is analytic in the complex variable s and only the 
second term can create poles at the complex values s — Sj where the condition 
( p^ is satisfied. We suppose here that these poles are simple. Near the pole 
s = Sj, we find a divergence of the form 



I - Q(s) S - S, . 3^Q(^^. 



^3 

Because of the definition (O), we have that 



X] • 9.Q(s, ).x,=--J2lb [b] * X3 [b] (43) 



V 

b 



In this way, we can identify the relaxation rates of the asymptotic time evolution 
of the physical averages with the roots of the characteristic determinant (p^. 
We can also identify the right eigenstates as 



(^1*.) - E ^^[^l f / ' x] dx (44) 

which is expected by the previous expression (|2^) for the right eigenstates, and 
the left eigenstates as 

From the definition (^6|) of the inner product, we infer that the right eigenstate 
associated with the resonance Sj is given by the following vector of 2B functions 

^,[b,x]^xAb]e-'^^ (46) 
while the corresponding left eigenstate is given by 

^ , \b x] = ^-^^^^ e*^* ^ (47) 

^^'^ Eb'k'Ub'hAb']* ^ ' 

which ends the construction of the spectral decomposition under the assumption 
that all the complex singularities of the Laplace transform of the Frobenius- 
Perron operator are isolated simple poles. 
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4.2 The classical zeta function 



The relaxation of the probabihty density is thus controlled by the relaxation 
modes which are given by the eigenvalues and the eigenstates of the Frobenius- 
Perron operator. As we said, the eigenvalues of the Frobenius-Perron operator 
are determined by the solutions {sj} of the characteristic determinant [see Eq. 




det [I - Q(s)] = 

These solutions are complex numbers which are known as the PoUicott-Ruelle 
resonances if they are isolated roots. 

We will rewrite Eq. (|lj) in a way that is reminiscent of the Selberg-Smale 
zeta function. With this purpose, we first consider the identity 

det(l - Q) = cxp tr ln(l - Q) 

and the expansion 

ln(l-Q) = -^^. 

n—l 

from which we infer the following identity 
det(l - Q) = exp 

Now, we have that 

(Q")66 = ^ Qbb„-i ■ ■ -Qb-ibiQbib 

6lfc2---6„-l 

and from Eq. (|l3|) we find 

trQ"=^(Q")bb- E Pbb^^,■■■Pb.MPb,be-^'^''r.-^+■■■+l^^+h) 

b bbib2...bn-i 

Note that this is a sum over closed trajectories composed of n lengths in the 
graph. The factor = Pbb„-i ■ ■ ■ PbibiPb^b plays the role of the stability factor 
of the closed trajectory 6&162 . . . 6n-i and following this analogy we define the 
Lyapunov exponent Ap per unit time as 

Al - CXp(-Aj,rp(")) = Pfc6„_, . . . Pb,b^Pb^b (48) 

where Tp"'' is the temporal period of this closed trajectory. We shall con- 
sider primitive (or prime) periodic orbits and their repetitions. A periodic 
orbit composed of n lengths can be the repetition of a primitive periodic or- 
bit composed of Up bonds if rt = rup and r is an integer called the repeti- 
tion number. With this definition the total period of the orbit is given by 



E 

n=l 



tr Q" 
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Tp^^ = {Ibi + ^2 + ■ ■ ■ + 'b„-i + h)/v = rlp/v with Ip the length of the primitive 
periodic orbit. Accordingly, we have the following relation for the Lyapunov 
exponent of the prime periodic orbit p — bib2 ■ ■ ■ bn^ composed of Up bonds: 



We can thus write 



trQ"= ^ npe-^-^^p e-^^^p 



(49) 



where we have explicitly considered the degeneracy Up of the orbit due to the 
number of points (vertices) from where the orbit can start. Accordingly, we get 



tr Q" 



E 

peP„ 



1 H 



The sum over n of is equivalent to the sum over all the periodic orbits and 
their repetitions, i.e.. 



n— 1 p r— 1 

We recognize in the right-hand side the expansion of the logarithm, whereupon 



In 



1 _e-(^''+") " 



We have thus 



det(l - Q) = cxp X! ^ ° 



(Ap + s)- 



or 



det(l - Q) = Hp 



1 _c-(^''+") " 



(50) 



which is the Selberg-Smale zeta function for the classical dynamics on graphs. 



4.3 The Pollicott-Ruelle resonances 

The zeros of the zeta function ( pO| ) are the so-called Pollicott-Ruelle resonances. 
The results here above show that the spectrum of the Pollicott-Ruelle resonances 
controls the asymptotic time evolution and the relaxation properties of the 
dynamics on the graphs. In general, the zeros Sj are located in the half- plane 
Resj < because the density does not grow exponentially in time. 

The spectrum of the zeros of the Selberg-Smale zeta function allows us to 
understand the main features of the classical Liouvillian time evolution of a 
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system. Let us compare the classical zeta function (^0|) for graphs with similar 
classical zeta functions previously derived for deterministic dynamical systems 



[ [12| 13 1 . For Hamiltonian systems with two degrees of freedom, the classical 
zeta function is given by two products: (1) the product over the periodic orbits 
as in the case ( |50| ) of graphs; (2) an extra product over an integer m = 1, 2, 3, ... 
associated with the unstable direction transverse to the direction of the orbit. 
This integer appears as an exponent of the factor associated with each periodic 



orbit [13|. As a consequence of this extra product, some zeros of the zeta 
function are always degenerate for a reason which is intrinsic to the Hamiltonian 
dynamics of a system with two or more degrees of freedom. Accordingly, Jordan- 
block structures are possible in typical Hamiltonian systems. 

In contrast, no such degeneracy of dynamical origin appears in classical 
graphs because no integer exponent affects the periodic-orbit factors in Eq. 
(|o|). In general, this property does not exclude the possibility of degenerate 
zeros which may appear for reasons of geometrical symmetry of a graph or for 
a particular choice of the parameter values defining a graph. However, such 
degeneracies are not expected for typical values of the parameters which are the 
transition probabilities Pbb' and the lengths Examples will be given below 
that illustrate this point. According to this observation, Jordan-block structures 
should not be expected in typical graphs. 

Different behaviors are expected depending on whether the graph is finite or 
infinite. 

4.3.1 Finite graphs 

Finite graphs are composed of a finite number of finite bonds. In this case, 
the matrix Q(s) is finite of size 2B x 2B with exponentials exp{—slb/v) in 
each element. The characteristic determinant and therefore the zeta function 
( pO| ) is thus given by a finite sum of terms with exponential functions of s. As a 
consequence, the Selberg-Smale zeta function is an entire function of exponential 
type in the complex variable s: 



\Z{s)\ < K exp[ — |s| 



where K is a positive constant and itot — Y^b=i is the total length of the 
directed graph (which is finite by assumption). Hence, the zeta function is 
analytic and has neither pole nor other singularities. In general, such a zeta 
function only has infinitely many zeros distributed in the complex plane s. 

The finite graphs form closed systems in which the particle always remains at 
finite distance without escaping to infinity. For closed systems, we should expect 
that there exist equilibrium states defined by some invariant measures. Such 
equilibrium states are reached after all the transient behaviors have disappeared 
in the limit t +oo. 



According to the spectral decomposition (37), the equilibrium states should 



thus correspond to vanishing relaxation rates sj — 0. Whether the equilibrium 
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state is unique or not is an important question. In the affirmative, the system is 
ergodic otherwise it is nonergodic. Because of the definition (|l3|), we have that 
Qbb' (0) = Pbbi so that the value s = is a root of the characteristic determinant 
( p^ if the matrix P of the transition probabihties admits the unit value as 
eigenvalue. Because of the condition we know that the unit value is always 
an eigenvalue of P. The question is whether this eigenvalue is simple or not. If 
it is simple, the equilibrium state is unique and the system ergodic otherwise it 
is multiple and the system nonergodic. 

In order to answer the question of ergodicity, let us introduce the following 
definition: 

P is irreducible m 'ib,b' , 3n : (P")fcb' > 

Then, we have the result that: 

The classical dynamics on a finite graph is ergodic if the matrix of the tran- 
sition probabilities is irreducible. 

Indeed, if the transition matrix is irreducible all the bonds are interconnected 
so that there always exist n transitions that will bring the particle from any 
bond 6' to any other bond b. It means that the graph is made of one piece, 
i.e., the dynamics on the graph is said to be transitive. The irreducibility of 
the transition matrix implies the unicity of the equilibrium state because of the 
Frobenius-Perron theorem |Q : 

// a matrix has non-negative elements and is irreducible, there is a non- 
negative and simple eigenvalue which is greater than or equal to the absolute 
values of all the other eigenvalues. The corresponding eigenvector and its adjoint 
have strictly positive components. 

We notice that the transition matrix P is non-negative and that no eigenvalue 
is greater than one because all the matrix elements obey < Pbb' < 1 and, 
moreover, Eq. (jj) holds. On the other hand, we know that the unit value is an 
eigenvalue also because of (0). Therefore, if the transition matrix is assumed to 
be irreducible, the eigenvalue 1 is simple. According to Eq. (|2^), the equilibrium 
state of relaxation rate Sq = is given by the unique positive eigenvector Xo = 
P ■ Xo = Q(0) ■ Xo corresponding to the simple eigenvalue 1 as 

*o[6, x] = xoib] , for < x < Zfc (51) 

The corresponding adjoint eigenvector of P is Xo[^] = Ij V6. The positive com- 
ponent Xo[^] of the right eigenvector gives the probability to find the particle 
in the bond b at equilibrium. These components obey the probability normal- 
ization X]fcXo[^] = 1- This equilibrium state defines an invariant probability 
measure in the space of trajectories: 

■ • ■ bibo) = Pb„_ifc„_2 • ■ ■ Pb2biPbiboXo[bo] (52) 

Examples of nonergodic graphs are disconnected graphs. 
The classical dynamics on a closed graph is said to be mixing if there is no 
pure oscillation in the asymptotic time evolution, i.e., if there is no resonance 
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with Re Sj — except the simple resonance sg = 0. According to Eq. we 
have for a mixing graph that 

lim {A)t - V xo[b] ^ A[b, x] dx (53) 

An example of a graph which is ergodic but non-mixing is a single bond of 
length g between two vertices. Its zeta fmiction is Z{s) = 1 — exp{—2sg/v) so 
that its resonances are 

Sj — i V — j with 7 G Z 
9 

Except So = 0, all the other resonances are pure imaginary so that the dynamics 
is oscillatory as expected. 

4.3.2 Infinite grapiis of scattering type 

Graphs of scattering type can be constructed by attaching semi-infinite leads 
c to a finite graph. These semi-infinite leads are bonds of infinite length. As 
soon as the particle exits the finite part of the graph by one of these leads it 
escape in free motion toward infinity, which is expressed by the vanishing of the 
following probabilities between the semi-infinite leads c and every bond b of the 
finite part of the graph: 

Ptc = V6 

and 

Pcb - Vfe 
Therefore, Qbc — Qtb — 0, V6, and 

det(l - Q) = det(l - Q) 

where Q is the matrix for the finite part of the graph without the scattering 
leads. Since the leads cause the particle to escape to infinity, the probability for 
the particle to stay inside the graph is expected to decay. Therefore, the zeros 
of the Selberg-Smale zeta function are located in the half-plane Re Sj < and 
there is a gap empty of resonances below the axis Res = 0: Resj < sq < 0. 
The resonance so with the largest (or smallest in absolute value) real part is real 
because the classical zeta function is real. This leading resonance determines 
the exponential decay after long times which we call the classical escape rate 
7(.i = — Sq (or the inverse of the classical lifetime of a particle initially trapped 
in the scattering region Td = l/7ci). The trajectories which remain trapped 
form what we shall call a repeller because it is the analogue of the repeller in 
deterministic dynamical systems with escape p^ . 

An invariant measure can be defined on this repeller by applying the Frobenius- 
Perron operator to the non-negative matrix Q(so) evaluated at the leading reso- 
nance. This matrix has a leading eigenvalue equal to one and the corresponding 
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left and right eigenvectors are positive. A matrix of transition probabilities on 
the repeller can be defined by 

Tlbb' = Xo[b] Qbw ^^TT (54) 

which leaves invariant the probabilities 

7r[6] = MlMl (55) 
Xo -Xo 

of finding the particle of each bond h in its motion on the repeller. These 
probabilities obey 

^Hfefc, = 1 , and ^Hfab, 7r[6'] = 7r[6] (56) 

b V 

and the invariant measure on the repeller is defined as 

■ ■ ■ ^i^o) — n6„_ib„_2 • • ■ nb2bjnf,jf,o7r[6o] 

As an example, consider the graph formed by one bond of length g that 
joins two vertices and two scattering leads attached to one of these vertices. 
The repeller consists here only of one unstable periodic orbit. Thus we look for 
the complex solutions of 



1 _e-(^p+^)v ^ 



that is 



27r 

Ap + s = — « w — i , with i G Z 

where 1^ = 2g and (IpXp/v) = — ln(l/9) which follows from Eqs. (|l|) and (p8|). 
Accordingly, we get 

ln9 TT , ^ 

Si = —V \- i V — j , with 7 e Z 

Therefore, all the resonances have the lifetime Tfi = = —f-^. This lifetime 
coincides with the quantum lifetime obtained from the resonances of the same 
graph Another system having this peculiarity is the two-disk scatterer 

[ p^ . This property is due to the fact that there is only one periodic orbit. In 
the presence of chaos and thus infinitely many periodic orbits, the quantum 
lifetimes are longer than the classical ones [nil n5[ . 



5 The chaotic properties 
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5.1 Correspondence with deterministic chaotic maps 

The previous results show that the classical dynamics on a graph is random. 
It turns out that this dynamical randomness is not higher than the dynamical 
randomness of a deterministic chaotic system. 

In order to demonstrate this result, we shall establish the correspondence 
between the random classical dynamics on a graph and a suspended flow on a 
deterministic one-dimensional map of a real interval. As aforementioned, the 
trajectories of the random dynamics on a graph are in one-to-one correspondence 
with bi-infinite sequences giving the directed bonds successively visited by the 
particle, • • • b-2b-ibobib2 • ■ ■ , which is composed of integers 1 < < 2B. For 
simplicity, we shall only consider the future time evolution given by the infinite 
sequence 60&1&2 • • ■ • With each infinite sequence, we can associate a real number 
in the interval < y < 2B thanks to the formula for the 2i?-adic expansion 

n=0 ^ ' 

Accordingly, the directed bond b' is assigned to the subinterval 6' — 1 < y < b' . 
Each of these subintervals is subdivided into 2B smaller subintervals 

Yb-i^b' <y < Yb,b' , 

with Yb,b' = n-i,b' + Pbb' , >o,b' = b' -I , and Fas.b' = b' (58) 

The one-dimensional map is then defined on each of these small subintervals by 
the following piece-wise linear function 

yn+l = (f){yn) = 77— [yn - 5^6-1. b') + 6 - 1 , for < yn < Yb^b' 

rbb' 

(59) 

Since the transition probabilities are smaller than one, < Pbb' < 1, the slope 
of the map is greater than one: 1 < As a consequence, the map ( |59| ) is in 
general expanding and sustains chaotic behavior. 

The suspended flow is defined over this one-dimensional map with the follow- 
ing return-time function giving the successive times tn of return in the surface 
of section: 

tn+l = tn + T(2/„) (60) 

with 

T{yn) = — , for Yb-i.b' < yn < n,6' (61) 

V 

For finite and closed graphs, the invariant measure of the one-dimensional 
map ( ^ is equal to 

Pcq(y) = Pb' , for b'-l<y< b' (62) 



23 



For infinite graphs of scattering type, the function (|9|) maps the subinter- 
vals associated with the semi-infinite leads outside the interval < y < 2B, 
generating an escape process. For such open graphs, the one-dimensional map 
selects a set of initial conditions of trajectories which are trapped forever in the 
interval < y < 2B. This set of zero Lebesgue measure is composed of unstable 
trajectories and is called the repeller. Typically, this repeller is a fractal set. 

We notice that, for closed graphs, an isomorphism can even be established 
between the dynamics in the space of bi-infinite sequences and a two-dimensional 
area-preserving map according to a construction explained elsewhere Q . 

5.2 Characterization of the chaotic properties 

The chaotic properties can be characterized by quantities such as the topological 
entropy, the Kolmogorov-Sinai entropy, the mean Lyapunov exponent, or the 
fractal dimensions in the case of open systems. All these quantities can be 
derived from the so-called "topological pressure" P{(3) (l^. This pressure can 
be defined per unit time or equivalently per unit length since the particle moves 
with constant velocity v on the graph. 

The topological pressure can be defined in analogy with the definition for 
time-continuous systems. For this goal, we notice that time is related to length 
hy V = l/t and that the role of the stretching factors is played by the inverses of 
the transition probabilities in the context of graphs. Accordingly, the topological 
pressure per unit time is defined by 

PiP) ^ hm J In {Pk-iK-. ■ ■ ■ Pbibof (63) 

ba---bn-2 
L<hg-\ hi6„_2<i+AL 

where the sum is restricted to all the trajectories that remain in the graph 
and do not escape (i.e., on the repeller) and that have a length that satisfies 

L < lbo-\ h lb„_2 < L + AL (cf. Rcfs. 0). The dependence on AL 

disappears in the limit L oo. 

The equation (|3|) can be expressed by the condition that the pressure is 
given by requiring that the following sum is approximately equal to one in the 
limit n ~f oo 

I'-n^oo E {Pb.-ib„-.---Pbibof e-i^WK-^+-+K) (64) 

bo---bri-2 

which is equivalent to requiring that the matrix Q(s; /3) composed of the elements 

Qtb'{s;P) = {Pwf e-''-^ (65) 

with s — P{0) has the eigenvalue 1 as its largest eigenvalue. As a consequence, 
the topological pressure can be obtained as the leading zero of the following zeta 
function 

Z(s;/3) =det[l-Q(s;/3)] (66) 
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or, equivalently, as the leading pole s = P{l3) of the Ruelle zeta function 

The different characteristic quantities are then determined in terms of the 
topological pressure function as follows : 

• The escape rate is given by 7ci = — P(l); 

• The mean Lyapunov exponent by A = — P'(l); 

• The Kolmogorov-Sinai entropy is determined by ft-KS = A — 7ci = P(l) — 
P'(l); 

• The topological entropy by /itop = ^(0); 

• The Hausdorff partial dimension of the repeller of the corresponding one- 
dimensional map (^9|) is the zero of P(/3), i.e., P(c?h) — 0. 

The mean Lyapunov exponent, the escape rate, and the entropies are defined 
per unit time. The mean Lyapunov exponent characterizes the dynamical insta- 
bility due to the branching of the trajectories on the graph. On the other hand, 
the KS entropy characterizes the global dynamical randomness. Both would be 
equal if the graph was closed and the escape rate vanished. We shall say that 
the dynamics on a graph is chaotic if its KS entropy is positive, h^s > 0. We 
emphasize that a dynamics with a positive Lyapunov exponent is not necessarily 
chaotic. A counterexample to this supposition is given by the open graph at 
the end of the previous subsection. The repeller of this graph is composed of 
a single periodic orbit and its Lyapunov exponent is equal to the escape rate: 
A = 7ci = (vln3)/(7. Accordingly, its KS entropy vanishes in agreement with 
the periodicity of this dynamics. 

We notice that the escape rate is related to the leading PoUicott-Ruelle 
resonance by 7ci — —sq. Indeed, when (3 — 1 the zeta function ( |66|) reduces to 
the previous one given by Eq. ( po|) which has the PoUicott- Ruelle resonances as 
its zeros. 

Moreover, we have the following properties: 

1. The topological entropy is independent of the transition probabilities Phb' 
of the graph. 

2. The Hausdorff dimension is independent of the lengths If, of the graph. 



The first property is deduced from Eq. ( pTj ) when we set /3 = to calculate 
the topological entropy. In this case, we observe that the Lyapunov exponents 
disappear from the zeta function ( |67| ) which thus depends only on the lengths 
of the periodic orbits of the graph. As a consequence, the topological entropy, 
which is given by the leading pole of the Ruelle zeta function ( |67| ) with (3 = 0, 
depends only on the lengths of the bonds. 
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Figure 1: The topological pressure for a fully connected pentagon with L — 10 
and L — 7Q leads attached to each vertex. The velocity is w = 1 . 



The second property can be inferred from Eq. ( pi] ) or, similarly, from the 
characteristic determinant (|6^) for the matrix (|65|). Indeed, since the Hausdorff 
dimension is the zero of the pressure function the lengths now disappear when 
we set s = P{f3) = in either Eq. ( |64[ ) or Eq. (|65|). Accordingly, the Hausdorff 
dimension depends only on the transition probabilities. 



6 Scattering on open graphs 

We shall consider some examples that illustrate the previous concepts. Consider 
the fully connected pentagon with L scattering leads attached to each vertex. 
Since the topological entropy /itop is independent of the Lyapunov exponents it 
is independent of the number of scattering leads attached to each vertex. This 
is observed in Fig. || where we depict the topological pressure for the fully 
connected pentagon. 

Moreover, we observe that the escape rate 7ci = — ^'(1) for the pentagon 
with L = 70 is smaller than the escape rate for the pentagon with L = 10. 

2 

This behavior has a simple interpretation. Since we use Phv = — <5p^, , 
the transmission probability from bond to bond decreases and the reflection 
probability increases as the valence of the vertex increases. Therefore, 
as the number of leads increases, a particle on the pentagon has a smaller 
probability to escape and a larger probability to be reflected back to the same 
bond. Accordingly, the escape rate diminishes. 

The examples of Fig. |^ shows that, indeed, the Hausdorff dimension is 
independent of the bond lengths. 
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Figure 2: The topological pressure for a fully eonneeted pentagon with L = 10. 
The curve with the crosses is obtained when all bonds as unit length and the 
curve with the circles for a set of incommensurate lengths. The velocity is t; = 1. 



As we see in these examples, the dynamics on typical graphs is characterized 
by a positive KS entropy /iks > 0- In this sense, the classical dynamics on typical 
graphs is chaotic. 



7 Diffusion on infinite periodic graphs 
7.1 The hydrodynamic modes of diffusion 

If the evolution of the density in an infinite periodic graph corresponds to a dif- 
fusion process, then the phenomenological diffusion equation should be satisfied 
in some limit. For instance, if the periodic graph forms a chain extending from 
X = — oo to a; = oo then on a large scale (much larger than the period of the 
system) the density profile should evolve according to the diffusion equation 



Let us notice that a; is a one-dimensional coordinate of position along the graph 
which is a priori different from the position along each bond. 
This equation admits solutions of the form 

Pk = exp [s{k) t] exp{ikx) 

with the dispersion relation 

s{k) = -Dk^ (69) 
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that relates the eigenvalue Sk to the wave number k. These solutions are called 
the hydrodynamic modes of diffusion. The inverse of the wave number gives the 
wavelength £ = 271 /k of the spatial inhomogeneities of concentration of particles. 

For a system such as a graph, we expect deviations with respect to the 
diffusion equation which only gives the large-scale behavior of the probability 
density and not the behavior on the scale of the bonds. Moreover, we may 
also expect the existence of other kinetic modes of faster relaxation than the 
leading diffusive hydrodynamic mode. In order to obtain a full description of the 
relaxation, we have to compute the eigenvalues of the evolution operator for an 
infinite periodic graph. One of these eigenvalues will have the dependence of Eq. 
( |69| ) for small k which allows us to obtain from it the diffusion coefficient of the 
chain. We shall start by considering periodic graphs in a d-dimensional space to 
show the generality of the method but then wc will specialize to periodic graphs 
that form one-dimensional chains. 

7.2 Fourier decomposition of the Probenius-Perron oper- 
ator 

In spatially extended systems which form a periodic lattice, spatial Fourier 
transforms are needed in order to reduce the dynamics to an elementary cell 
of the lattice In this reduction a wave number k is introduced for each 
hydrodynamic mode. The wave number characterizes a spatial quasiperiodicity 
of the probability density with respect to the lattice periodicity. Indeed, the 
wavelength of the mode does not need to be commensurate with the size of a unit 
cell of the lattice. Each Fourier component of the density evolves independently 
with an evolution operator that depends on k. Accordingly, the PoUicott-Ruelle 
resonances will also depend on k. We shall implement this reduction starting 
from the master equation ( p5| ) of the fully periodic graph and construct from it 
the evolution operator for the unit cell. A few new definitions are needed before 
proceeding with this construction. 

The periodic graph is obtained by successive repetitions of a unit cell. Such 
graphs form a Bravais lattice C. A lattice vector am is centered in each cell 
of the Bravais lattice with m = (mi,m2, ...,md) E Z"*. The lattice vectors are 
given by linear combinations of the basic vectors of the lattice 

am = "Maioo---oo + '7i2aoio...oo H ^ ?"daooo---oi e C 

(for a one-dimensional chain we have d = 1 and am — vri eTL). We shall split the 
coordinate [6, x] which refers to an arbitrary bond in the infinite graph into the 
new coordinates ([b, a;],a) where the first pair refers to the equivalent position 
in the elementary unit cell to which the dynamics is reduced. That is, the 
bond b is associated with h and x is the position in that bond. The third term 
represents a vector in the Bravais lattice that gives the true position of the bond 
h with respect to the position of the original unit cell, that is h is obtained by 
translating the bond b to the cell in the Bravais lattice identified by the vector 
a. We may introduce the notation h = Talb) where the translation operators 
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Ta assign to a bond b the corresponding bond in the unit cell characterized by 
the vector a. 

Accordingly, the density in the graph is represented by a new function |^ p 
related to the old one by /9([b, a;], a, t) — p{\b, x\,t). 
We define a projection operator by 

i^k = ^ exp(-ik • a)S"', (70) 

in terms of the spatial translation operators 

S"'/([b,a;],a',i) = /([b, x], a' + a, t) for a, a' e C. 

The projection operator (|7^) involves the so-called wave number k. This later 
is defined on the Brillouin zone B of the reciprocal lattice £. The volume of the 
Brillouin zone is 



|det(aioo---oo, aoio---oo, aooo---oi)| 
The operators ( |70| ) are projection operators since 

E^E^, ^ \B\ <5(k - k')^k 
which is a consequence of the relation 

|i|^exp(^k.a)= ^ 5(k - k'). 

The identity operator is recovered by integrating the projection operator over 
the wave number 

If p is the density defined on the infinite phase space, the function i?kP is 
quasiperiodic on the lattice 

i?kp([b, x], a, t) = ^ exp(— ik • a')/9([b, x], a + a', t) 

a'G£ 

= exp(ik • a) exp(— ik • a )p([b, a;],a , i) 

a"G£ 

= exp(ik • a)£'kp([b, x], 0, t) — exp(ik • a)/3k([b, x\,t). 

We have therefore a decomposition of the density over the infinite phase space 
into components defined in the reduced phase space and which depends contin- 
uously on the wave number k. 

^ We keep calling p this new density. 
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Consider the master equation ( p5| ) of the full periodic graph. Applying at 
both sides the operator wc get from the quasiperiodicity of E-^^p [see the 
previous equation] 

exp(ik ■ am)pk([b,a;], i) = ^ Pw exp(ik • am')Pk ( [b',x], i 

b' \ ^ 

(71) 

where h = Ta„(b) and b' = Ta_^, (b') and also lb = Ih- Now, the translational 
symmetry implies that 

Pbb' ^ PT^^{h),T^ ,(b')=^Ta _ ,(b),b' 



Pk([b,a;],t) = X!X!-^^-„_„'(b),b' exp(-ik- am-m') Pk ([b',x],i 



Thus Eq. (|71|) is 

b' m' 

The time evolution of the k-component is therefore controlled by the matrix of 
elements 

Pbb'(k) = ^7'e,„(b),b' exp(-ik • a„,) (72) 

m 

We observe that, for a graph, there is at most one term in the sum of Eq. ([72|). 
Indeed, the coefficient ^^^^(bj.b' does not vanish if and only if the bonds ^^^^(b) 
and b' are connected with the same vertex of the infinite graph. Furthermore, 
there is one and only one translation for which Ta^ (b) and b' are connected 
with the same vertex. Accordingly, there exists a unique lattice vector a(b, b') 
such that 

Pbb'(k) = PT,„.„,(b),b' exp [-ik ■ a(b, b')] (73) 

Thanks to this matrix, we have that each Fourier component of the density 
evolves with an equation 

Pk([b,x],i) = 5j]Pbv(k) pu (^[b',x],i- ^±i^l^^ (74) 

Here the sum is carried out over all the directed bonds of the unit cell and the 
matrix P(k) of elements Pbb'(k) defined by Eq. ( [73| ) is a square matrix with 
dimension equal to the number of directed bonds in the unit cell. 

7.3 The eigenvalue problem and the diffusion coefficient 

To study the eigenvalues of t he F robcnius-Pcrron operator i?^ defined by Eq. 
( [74I ) we proceed as in Subsec. \i.2\ We introduce the following definitions 

t)dt 



PKs[b]^ e'^Vk([b,0], 
Jo 
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and 



Qbb'(s,k) = Pbb'(k) exp (-s^) (75) 

for the elements of the matrix Q(s,k). Then, taking the Laplace transform of 
Eq. we get 



Pk,s [b] = ^(5bb'(s,k) /9k,s [b'] 

b' 



which has a solution only if the following classical zeta function vanishes 

Z{s; k) = det[l - Q(s, k)] = H " ^ q (yg) 

p 

where the product extends over the prime periodic orbits p of the unit cell of 
the graph and where ap = X]r=i ^(bz+ijbi) is the displacement on the lattice 
along the periodic orbit p = bib2 • • • hn^ of prime period rip.^ From Eq. ([76|), we 
obtain the functions Sj(k) and the corresponding eigenstates XjM which satisfy 

X3M = Qh(k),k] -XjM- 

Eq. ( [76| ) shows that the PoUicott-Ruelle resonances sj (k) are the zeros of a new 
classical Selberg-Smale zeta function defined for the spatially periodic graphs. 
The eigenvalues and eigenstates of the Frobenius-Perron operator i?^ defined 



by Eq. (|7j) are constructed in the same way as in Subsec. |3.2| . If we denote by 
\l/j_k[b, 2:] the eigenstates of the flow, the results are 



with 



i?^*,,k[b,x] =e^^('^)*^',-k[b,a;] 



*j,k[b, x]=e Xj,k[b] for < a; < Zb 



where Sj(k) is a solution of Eq. (|76|). 

For k = 0, Eq. ( [74[ ) represents the evolution of the density in the unit cell 
with periodic boundary conditions. The periodic boundary conditions transform 
two bonds in a loop. In this way, for k = 0, we are studying the evolution of 
a closed graph and we saw in Subsec. |3.2| that, for a closed graph, the value 
s = is a solution of Eq. ( [76| ) and the associated eigenstate corresponds to the 
invariant measure or equilibrium probability which is given by the 2i?-vector 
{B being the number of bonds in the unit cell) 

Xo = ^(l,l,... ,1). (77) 



''Here, a(bi+i,bi) denotes the jumps over the lattice during the transition between the 
bond hi and the bond bi_|_i. We have to consider that a(bi_|_i,bi) = for transitions between 
bonds in the same unit cell. 
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We may thus expect that, for k small enough, there exists a zero so(k) of Eq. 
7m and a corresponding eigenstate 



Xo,k = Q [so(k),k] -xo.k 



(78) 



such that so(k)g^;^0. This particular resonance can be identified with the dis- 
persion relation of the hydrodynamic mode of diffusion since this latter is known 
to vanish at k = as 



SO that the diffusion matrix is obtained as 

n 1 5'«o(k) 

Uat3 - - 77 



2 dkadkp 



(79) 



k=0 



7.4 Example of infinite graph with diff"usion 

We illustrate the results of this section with a very simple example of a one- 
dimensional chain. Consider the graph used in Subsubsec. 4.3.2. This graph of 
scattering type can be used as a unit cell for a periodic graph. The right lead 
is connected with the left lead of an equivalent graph and so on. This graph 
looks like an infinite comb. The unit cell can be considered as composed by two 
bonds, say b and a. The bond b connects the dead vertex 2 with the vertex 1 
and the bond a connects the vertex 1 with the vertex 1 of the next cell. Thus 
the valence of the vertices are V2 — 1 and vi = 3 respectively. The transition 
probabilities Pbb> are given by 



bb' 



\Tbb' 



with, P,f, = l,P,t, = |,R 



bb 



Pi. 



and zero otherwise. 



We first construct the matrix P(fc). We note from its definition in Eq. (73) that 
^hh'{k) = Pbb' for bonds that belong to the unit cell. The fc-dependent factors 
come from bonds that connect consecutive cells. These bonds are: 



The bond a of the cell at the left-hand side characterized by am = — 1 is 
connected with the bonds a and b. This gives the contributions Paa(A) = 
|e+^'= and P^Jfc) = |e+'^ 

The bonds a and b of the cell at the right-hand side characterized by 
am, = +1 are connected with the bond a. This gives the contributions 



Paa(fc) = le-'" and Pab(fc) = fe- 
Therefore, P{k) is the 4x4 matrix with entries 



P{k) = 



-ik 







9 

4p-ifc 
9*^ 
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Figure 3: The first two branches s{k) of Poliicott-Ruelie resonances of the infi- 
nite comb graph, as obtained from Eq. (^o|). The branch containing the origin 
s — and fc = is the dispersion relation of the diffusive mode. The other 
branch is associated with a kinetic mode of faster relaxation. 

where the columns and rows are arranged in the following order (b, b, a, a). The 
matrix Q(s, k) is obtained by multiplication with the diagonal matrix (5bb'e~ . 
The determinant in Eq. ( j?^ ) can be computed and gives 



where a is the length of the bond a and g is the length of the bond b. As we 
said the solutions of det [I — Q(s, k)] = gives the desired functions Sj{k). For 
this example, we plot the first branches in Fig. ^ where we observe that, indeed, 
only one branch includes the point s = at fc = 0. This unique branch can be 
identified with the dispersion relation of the hydrodynamic mode of diffusion. 

The diffusion coefficient is obtained from the second derivative of the first 
branch at fc = 0. This can be analytically computed for this particular example 
as follows. We consider s « Mm{v/a,v/g} and fc << 1 and expand Eq. (pO|). 
After some simple algebra, we get 



from which we obtain that the diffusion coefficient defined by Eq. (179) is 



det [l-Q(s,A;)] = l + |=e 




-) 

(80) 



cos k 




V 



(81) 



D = 



a + g 



33 



In general, the diffusion coefficient has the units of [I/^]/[T]. This is because 
the wave number has the units of Since we have considered as a di- 

mensionless number, the diffusion coefficient has here the units of In 
this example the standard units can be recovered by considering k = na with a 
standard wave number k where a is the period at which the unit cell is repeated. 
In these units, we would have obtained D = v-^. But, in general, for a more 
complicated graph, there is no bond length that we can associate with the peri- 
odicity of the chain as in this example, which is the reason why we consider the 
dimensionless parameter k. In this sense we are considering the space in units 
of the unit cell of the periodic chain. 

7.5 A Green-Kubo formula for the diffusion coefficient 

In the previous example, we have seen that the diffusion coefScient is inversely 
proportional to the total length of the unit cell. This is a general property that 
follows from a general expression for the diffusion coefficient that we shall now 
obtain. From now on we shall consider one-dimensional chains]^. Accordingly, 
fc is a scalar wavenumber and no longer a vector. 

Consider the vector Xo,k defined by Eq. and the eigenvector Xo.fc of the 
adjoint matrix Q [so(fc), fc]^ defined by 

Q[so(fc),fc]^Xo,fc=Xo,fe (82) 

or equivalently 

xU=xU-QN(fc),fc] (83) 

This eigenvector satisfies 

Xo,fc[b] 1 (84) 

Such vectors are normalized as 

(Xo,fc:Xo,fc) ^^XoAWxoAM = 1 (85) 

b 

On the other hand, due to Eqs. (|7^ ) and (^2|), we have 

(Xo,fe,Q(so,fc)-Xo.fe) -1- (86) 

with So — so{k)- Differentiating Eq. ( |85| ) and Eq. ( |8^ ) with respect to k we get 
respectively 

fdXo.k \ , f~ dxo.k\ „ 

The theory developed here is trivially extended for graphs that display periodicity in 
higher dimensions by considering the appropriate dimensionality for the vector a. 
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and 



dk 



iXO,k 



dk 



Q{so,k) 



■ Xo,fe 



dk 







where we have used Eqs. ( |78| ) and (|8: 
d 



dk 



Q(so,fc) 



Xo,fe ) = Xl^o.'^-W 

b,b' 



These last two equations imply 
d 



dk 



Qhh'iso, k) 



XoAb']=0 (87) 



Now we compute the derivative of Q. From Eq. (iTa) and Eq. ((73), we have 



^b' \ ^ . , N ^b' 

So — - Qhh'(so,k) — 

V J ~' 



Inserting this result into Eq. 



dsp 
dk 



exp 

|) and taking the limit k 
j_ dP 

2B ^ 



V dk 
0, we obtain 



^b,b' 



dk 



k=0 



k=0 



2B Sb,b' P 



bb' |fc=o ■ 



where we have used the fact that the limit k ^ Q implies So(fc) 0, and that 
Xo,fe[b'] ^ 2^'Xo,fe[b] ^ 1 because of Eqs. (^) and (||). Since Y.h Pbb'lfe^o = 
1, Vb', this reduces to: 



dk 



k=0 



b,b' 



dk 







(88) 



k=t) 



where Luc = X^b total length of the unit cell. The last equality (= 0) 

follows from the fact that the unit cell is connected with the neighboring cells 
in a symmetric way. For instance, in a one-dimensional chain, the "fluxes" from 
the left-hand side equal those from the right-hand side and the derivative with 
respect to k drops a sign that makes the sum vanishing. The reader can verify 
this property in the previous example of the comb graph. 

To obtain the diffusion coefficient we need the second derivative of so(fc). 
Therefore, we differentiate Eq. (|^ with respect to k and we evaluate at fc = 0. 
After some algebra and using Eq. (|8^), we get 



dfc2 



fe=0 



b.b' 



fe=0 



C^Pbb' 

b.b' 



dk 



2B 



dXQ,k\^'] , dxo,k\^V 



dk 



dk 



k=Q 



The explicit form for the diffusion coefficient is obtained from Eq. (^9|) if 
we compute the first derivatives of the eigenstates. We can write the equations 
that these quantities satisfy. In fact taking the derivative with respect to k of 
Eqs. (H) and (||) we have 



rfXo.fc[b] 
dk 



E 

b' 



dQhh' 
dk 



Xo,fc[b] 



E 

b' 



Qbb' 



dXo,k[^'\ 
dk 
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and 



= 2^Xo,fc[b] —Ti-+2^ Qbb' 



dk ^-^ ' dk ^-^ dk 

b b 

whose solutions are 

= 2. (I - Q)bb' -^xoMh 

b'b" 



dk 



b'b" 



dk 



In the Hmit k 0, these solutions can be written as 



dk 



k=0 



^ oo 

^ ^ 

b' b" n=0 



rfPb' 



b',b" 



b,b' 



dk 



fc=0 



and similarly 



dk 



EE 



dFu, 



b'.b" 



dk 



(P")b"b 



^==0 b',b"n=0 

Thus, the second term of Eq. becomes 



fe=0 



E 

b,b' 



dF 



bb' 



2J5 



c?Xo.fc[b'] , rfXo.k[b]* 



dk 



k=0 



-2 E E 

b,b',b",b"' n=0 



rfPbb' /p„s dPb"b"' 
ib'.b" 



dk 



dk 



k=0 



To evaluate and interpret this result we have to transform these expressions. 
First, we have to consider the derivatives of P(A;). This matrix is defined in Eq. 
([73|). Since only the nearest neighbors are connected, the lattice vector of the 
jumps can take only the values a(b, b') = 0, ±1 whether the particle crosses the 
boundary of the unit cell to the right-hand cell (+1), or the left-hand one (-1), 
or it stays in the same cell (0) during the transition b' — > b. Therefore 



Pbb'(fc) = PT„(,,,,,(b).b' e- 
The derivatives of this matrix are thus 

dPbb' 



-ika{h,h ) 



dk 



= -i a(b,b') P 



and 



~dk^ 



bb' 



-a(b,b' 



l\2 



bb' 



bb'- 
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Accordingly, the diffusion coefficient is given by 
1 d^so 



D 



2 dk^ 



k=0 



2L 



-y 



b.b' 



a(b,b')'P 



bb' 



2 E E«(b,b')Pbb'(P")v,b"a(b",b"')Pb"b" 



b"b"' n=0 



(90) 



fe=0 



In order to interpret this formula, we have to remember some definitions. If an 
observable is defined over M successive bonds, its mean value over the equilib- 
rium invariant measure of the random process is given by 



{A)= lim V A{bM-i---bibQ) ^l{bN■■■b^N) 

N^OO ^ — ^ 



(91) 



If the observable A depends only on two consecutive bonds as it is the case for 
the jump vector a(b,b'), its mean value takes the form 

(A) = J2 ^(^l^o) /^(^l^o) = J2 ^(^l^o) Pb.bo Xo[bo] = ^ E ^(^l^o) Pb^bo 



bobi 



bobi 



because of Eqs. ( |5^ ) and (77). According to the general definition (pi|), the 
time-discrete autocorrelation function of a two-bond observable is given by 

{A„-,Ao) = E A{b„i+ib„i)A{bibo) fi{b„i+ibm ■ ■ ■ bibo) 

bo---b,„ + i 

Because of Eqs. ( js^ ) and (^ again, we get 

{A^Ao) = ^ E Mbn^+lb„^)Pb^+,b^{P"'-'),^^,^A{b,bo)Pb,b, (92) 

bobibmbm + 1 

The terms of Eq. ( pO| ) are precisely of the form of Eq. ( |9^ ) with m = for 
the first term and m = n -I- 1 for the following ones, and with the observable 
A = a. Since the process is stationary we have that (a„ao) = (aoa_„) = (a_„ao) 
where the last equality follows from the commutativity of the quantities ao 
and a_„. Therefore, the term with the sum over n in Eq. ( pO[ ) is equal to 



2Bj2n=-oo risio i'^oa.n)- It is now clear that Eq. (90) for the diffusion coefficient 



IS 



n— — oo ^ ' n— — oo 

where (^) — Luc/{2B) is the mean bond-length of the unit cell and a„ = 0, ±1 
is the jump from one cell to another undergone by the particle in motion on the 
infinite graph. 
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Eq. (|9^) is nothing else than the Green-Kubo formula for the diffusion 
coefficient. If we define 



as the velocity along the x-axis that contributes to the transport, where Ai = 
(/) /v, we can write Eq. (|9^ ) in the more familiar Green-Kubo form 

\ I {vMvx{t))dt. 

In the time-discrete form ( |93| ) , we obtain the result that the diffusion coeffi- 
cient is proportional to the constant velocity v and inversely proportional to the 
mean bond-length of a unit cell. The diffusion coefficient is also proportional to 
the sum of the time-discrete autocorrelation of the jump a from cell to cell. 



8 Escape and diffusion on large open graphs 

In this section, we shall study the PoUicott-Ruelle resonances of open graphs 
characterized by a unit cell which is repeated a finite number of times. The 
particular example that we consider is depicted in Fig. IJ. We shall focus on the 
leading resonance that determines the escape rate from the system. We shall 
show that, for large enough chains (i.e., made of several unit cells), the classical 
lifetime corresponds to the time spent by a particle that undergoes a diffusive 
process in the chain before it escapes. 
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Figure 4: Open graph forming a chain composed of = 6 identical unit cells. 
The five lengths that compose the unit cell take different values. (In this regard, 
no specific length can be associated with the periodicity of the chain.) 

For the graph of Fig. ^ the transition probabilities from bond to bond Pf,f,/ 
are given by 

if the particle is reflected, i.e., 6 = 6'; 
^ for bonds h ^ b' which are connected ; 
otherwise . 

We have computed the spectrum of PoUicott-Ruelle resonances for different val- 
ues of the number N of unit cells. The leading resonance controls the asymptotic 
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Figure 5: Topological pressure for the chain of Fig. |4| but with N = 7. From 
this function, we get that Hks ~ 2.9 > and ft-top = 3.2330. 

decay. Since the leading resonance is isolated and at a finite distance from the 
real axis the decay is exponential as we explained, that is 

p{t) ~ exp (-7cii) 

where 7ci is the leading resonance, i.e., the escape rate. This is the generic 
behavior of the density in a classically chaotic open system and we refer to this 
as the classical decay. 

In Fig. H we depict the topological pressure for the chain of Fig. ^ with 

= 7 unit cells. 

When the chaotic dynamics is at the origin of a diffusion process the escape 
rate is inversely proportional to the square of the size of the system, more 
precisely the following relation is expected to hold 

7ci(iV) (94) 

with D the diffusion coefficient. ^ As explained in Sec. ^, the escape rate is 
related to the mean Lyapunov exponent and the KS entropy of the open chain 
of size N according to 

7ci(iV) = A(7V)-/iKsW (95) 

* This relation is obtained by solving the diffusion equation (tesl) in a system of size with 
absorbing boundary conditions at the borders, i.e, p(0) = and p(iV) = 0. The mode with the 
slowest decay is then given by sin(fcx) with k = and, from the dispersion relation (|6S|), we 
get Eq. ( ^). For large systems, when TV — ► oo, or equivalently fc — ► 0, D(N) = "/cii^)iN/iT)^ 
must approach the diffusion coefficient. 
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As a consequence of Eq. ( p^ , we find also for open graphs a known relationship 
between the diffusion coefficient and the chaotic properties 

D = lim — [X{N) - hKs{N)] (96) 

N^oo TT 

Here again D is in units of [i]/[T] because we did not associate a length with 
the period of the chain and thus the space is in units of unit cell. 





Figure 6: Lifetime of the chain of Fig. ^ as a function of the square of its 
size TV: (a) For iV = 3, 4, 5, 6, 7, 8, the slope of the linear regression is 0.247403 
and Eq. (|94| ) gives an approximate diffusion coefficient D ~ 0.4095. (b) For 
N = 46,56,66,76,86,106, the slope of the linear regression is 0.1946247 and 
Eq. (^) gives a better approximation D ~ 0.5205 for the diffusion coefficient 



We have computed the escape rate for chains of different sizes. If the dynam- 
ics indeed corresponds to a diffusion process, then Eq. (94 ) should be verified. In 
figure ^, we plot the classical lifetimes Td as a function of 7V^. Since Td = l/7ci) 
we observe the dependence of 7ci on N expected from Eq. (|93) . We may con- 
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elude from this result that the classical dynamics in the open chain is the one 
of a diffusion process. 
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Figure 7: The first branches s{k) of PoUicott-Ruelle resonances of the infinite 
graph corresponding to the open graph of Fig. ^. Here again we observe that 
there is only one branch that is identified with the hydrodynamic mode of dif- 
fusion, so{k), because it vanishes for fc = 0. 

The diffusion coefficient of the infinite graph can be obtained as we explained 
in Sec. ^ We depict in Fig. ^the leading and other PoUicott-Ruelle resonances 
of the infinite graph obtained by numerical calculation as a function of the 
dimensionless wave number k. The diffusion coefficient is given by the second 
derivative of the leading resonance so{k) evaluated at fc = 0. In this way, we 
obtain the numerical result 

D = 0.5318. (97) 

Accordingly, the diffusion coefficient of the infinite chain gives a reasonable 
estimate for the proportionality coefficient between 7ci and l/N"^ for the small 
chains of Fig. ^ and is a very good estimate for the large chains of Fig. 

In Fig. ^, we show how the effective diffusion coefficient D{N) = jci{N){N/7r)^ 
converges to the diffusion coefficient of the infinite chain D, as the chain becomes 
longer and longer {N oo). 

9 Diffusion in disordered graphs 

In a recent work, Schanz and Smilansky have considered the problem of 
Anderson localization in a one-dimensional graph composed of successive bonds 
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Figure 8: The effective diffusion coefficient D{N) = 7ci(iV)(iV/7r)^ as a function 
of N for the open graphs of Fig. ^. 

of random lengths with random transmission and reflection coefficients at the 
vertices. The classical dynamics corresponding to this quantum model defines a 
kind of Lorcntz lattice gases as studied in Refs. |2l[ . Indeed, these references 
describe Lorentz lattice gases consisting of a moving particle traveling with 
allowed velocities ±v on a one-dimensional lattice of scatterers. If the particle 
arrives at a scatterer it will be transmitted or reflected with probabilities p and 
q = 1 — p respectively. If the scatterers are randomly distributed the model 
describes the classical dynamics of the model by Schanz and Smilansky with 
identical transmission and reflection coefficients at all the scatterers. 

The classical dynamics of this model can be analyzed with the methods 
developed in the present paper, which provides the relationship with the cited 
works on the Lorcntz lattice gases. Using the methods of Sec. |[ we can write 
down the infinite matrix Q(s). The eigenstates of this matrix corresponding to 
the unit eigenvalue can be obtained by iteration along the chain according to 
x[b] = exp(s4/2w)wf, and x[b] — exp(sZf,/2w)uj| with 



(98) 



where b Z, and the matrix is defined by 

Mb 



r1/2 w r1/2 



(99) 
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with 



and 



96 

Pb 

_ _ 9t 
Pb ' Pb 



= £ 1 .'k (ICO) 



We notice that det = 1. If the chain was periodic If, ~ I, we would obtain the 
diffusion coefficient D = vlp/{2q) by assuming that Wb+i — exp{iKl)ui, in Eq. 
(p8[). In the dilute gas limit, the mean- field diffusion coefficient for the random 
graph is then given by replacing the bond length I by the mean bond length {I), 
leading to: 

^mf = v{l)^ (102) 
Zq 

For a disordered chain with N scatterers, the Pollicott-Ruelle resonances 
can be obtained by finding the resonances s for which the following equation is 
satisfied: 

N 



UN+l J ' V "1 

^ 0—1 ^ 



nM,(,s). (103) 



If the chain closes on itself, we must impose the periodic boundary conditions 
u-ffj~^ = Ui and utv+i = ui. If the chain is open and extended by two semi- 
infinite leads, we must consider the absorbing boundary conditions w^^^— ^ = 
and ui = 0. 



In Ref. 1 20 1, Ernst et al. have characterized the chaotic properties in such 
open graphs thanks to the escape-rate formalism by computing the topological 
pressure function of Sec. ||. In Ref. Appert et al. showed that the 

spatial disorder is at the origin of a dynamical phase transition associated with 
a singularity in the pressure function of the infinite disordered chain. 



10 Conclusions 

In this paper, we have introduced and studied the random classical dynamics of 
a particle moving in a graph. We shall show elsewhere ||ll|] that the dynamics 
here studied is the classical limit of the quantum dynamics introduced in Refs. 

We have shown that the relaxation rates of the time-continuous classical dy- 
namics can be obtained by a simple secular equation which includes the lengths 
of the bonds and the velocity of the particle. This secular equation has been di- 
rectly related to the eigenvalue problem of the time-continuous Frobenius-Perron 
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operator. The secular equation can be written as a classical zeta function de- 
fined as a product over the periodic orbits of the graphs. In this way, we have 
been able to define the relaxation rates, as well as chaotic properties such as the 
Lyapunov exponents and the entropies as quantities per unit of the continuous 
time. The chaotic properties are derived from a pressure function defined for 
each graph. 

For infinite periodic graphs, we have shown how to construct the hydrody- 
namic modes of diffusion and to compute a diffusion coefficient. Here also, the 
relaxation rates of the hydrodynamic modes are given by the zeros of a classical 
zeta function. Moreover, a Green-Kubo formula for the diffusion coefficient has 
been deduced from the eigenvalue problem for the Frobcnius-Pcrron operator of 
the classical dynamics on the graph. 

When the chain is open by considering a finite segment connected with scat- 
tering leads, the particle escapes after a diffusion process. In this case, we have 
computed the lifetime of the metastable states. This classical lifetime is given 
by the inverse of the escape rate which is related to the diffusion coefficient. Ac- 
cordingly, a known relationship between the diffusion coefficient and the chaotic 
properties jl^ is extended to the random classical dynamics on graphs. The 
case of infinite disordered graphs has also be considered. 

The interest of these results lies notably in the fact that the classical quan- 
tities here computed can be compared to equivalent quantities defined for the 
corresponding quantum problem, as shown elsewhere pl| . 
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